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Abstract 

We present an improved measurement of CP-violation parameters in B° — > (f>K s , K + K~ K s , 
and n'K s decays based on a 140 fb _1 data sample collected at the T(45) resonance with the Belle 
detector at the KEKB energy-asymmetric e + e~ collider. One neutral B meson is fully reconstructed 
in one of the specified decay channels, and the flavor of the accompanying B meson is identified 
from its decay products. CP- violation parameters for each of the three modes are obtained from 
the asymmetries in the distributions of the proper-time intervals between the two B decays. We 
find that the observed CP asymmetry in the B — > 4>K S decay differs from the standard model 
(SM) expectation by 3.5 standard deviations, while the other cases are consistent with the SM. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the standard model (SM), CP violation arises from an irreducible phase, the 
Kobayashi-Maskawa (KM) phase jj], in the weak-interaction quark- mixing matrix. In partic- 
ular, the SM predicts CP asymmetries in the time-dependent rates for B° and B° decays to a 
common CP eigenstate fcp\M- Recent measurements of the CP-violation parameter sin 20! 
by the Belle 0] and BaBar jj collaborations established CP violation in B° — > J/ipK s and 
related decay modes , which are governed by the b — > ccs transition, at a level consistent 
with KM expectations. 

Despite this success, many tests remain before one can conclude that the KM phase is 
the only source of CP violation. The B° —>■ <pK s decay, which is dominated by the b — > sss 
transition, is sensitive to new CP- violating phases from physics beyond the SM p. The 
other charmless decays B° —>■ K + K~K S and B° —>■ rj'K s , which are mediated by b — > sss, 
suu and sdd transitions, also provide additional information. Since the SM predicts that 
measurements of CP violation in these charmless modes should also yield sin20i to a good 
approximation 0,0], a significant deviation in the time-dependent CP asymmetry in these 
modes from what is observed in b — > ccs decays would be evidence for a new CP- violating 
phase. 

In the decay chain T(4S f ) — > B°B° — > fcpf tag, where one of the B mesons decays at 
time tcp to a final state fcp and the other decays at time t tag to a final state ft ag that 
distinguishes between B° and B°, the decay rate has a time dependence given by [2j 

e -\At\/T B J" n 

V(At) = {l + q- Ssm(Am d At) + Acos(Am d At) \, (1) 

where r^o is the B° lifetime, Am d is the mass difference between the two B° mass eigenstates, 
At = top — t tag , and the 6-flavor charge q = +1 (—1) when the tagging B meson is a B° 
(B°). S and A are CP- violation parameters; to a good approximation, the SM predicts 
S = sin 2^i, where £/ = +1(— 1) corresponds to CP-even (-odd) final states, and .4 = 
for both b — > ccs and b — > sss transitions. The present world-average S value obtained from 
previous measurements by the Belle [j| and BaBar 0| collaborations is within 2a of the 
SM expectation for B° rj'Kg and K + K~Kg, while a 2.7a deviation exists for B° — > (f)Kg. 
Measurements with a larger data sample are required to resolve this difference. 

Belle's previous measurement for B° — > <pKg, K + K~K$ and rj'Kg was based on a 78 
fb _1 data sample containing 85 million BB pairs. In this Letter, we report an improved 
measurement incorporating an additional 62 fb _1 for a total of 140 fb" 1 (152 million BB 



pairs). At the KEKB energy-asymmetric e + e~ (3.5 on 8.0 GeV) collider [111] , the T(4S') is 
produced with a Lorentz boost of (3 r j = 0.425 nearly along the electron beamline (z). Since 
the B° and B° mesons are approximately at rest in the T(45) center-of-mass system (cms), 
At can be determined from the displacement in z between the fcp and /tag decay vertices: 
At ~ (z C p - z tag )/ 'Plc^Az/ 'Prfc. 

The Belle detector [12] is a large-solid-angle spectrometer that includes a three-layer sil- 
icon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), time-of-flight (TOF) scintillation counters, and an 
electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect K\ mesons and to identify muons (KLM). 

We reconstruct B° decays to 4>K S and rj'K s final states for £j = —1, and B° —>■ K + K~K S 
decays that are a mixture of = +1 and —1. K + K~ pairs that are consistent with 
— > K + K~ decay are excluded from the B° — > K + K~K S sample. We find that the 
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K + K~Kg state is primarily = +1; a measurement of the = +1 fraction with a 140 fb _1 
data set yields 1.03 ±0.15(stat) ±0.05(syst), which is consistent with the previous result 
In the following determination of S and A, we fix £j = +1 for this mode. The intermediate 
meson states are reconstructed from the following decay chains: r/ — > p°(— > 7r + 7r~)7 or 77' — > 
7r + 7r _ ^(— > 77), iTg — > 7r + 7r~, and — > K + K~. Candidate Kg — > 7r + 7r~ and — > K + K~ 
decays are selected with the same criteria as those used for the previous branching fraction 
measurements 0|. The K% selection is slightly changed from the previously published CP 
asymmetry measurement |9j to improve the Kg purity. 

Since the meson selection is effective in reducing background events, we impose only 
minimal kaon-identification requirements. We use more stringent kaon-identification re- 
quirements to select non-resonant K + K~ candidates for the B° — » K + K~K S decay &}. We 
reject K + K~ pairs that are consistent with D° — > K + K~ , Xco ~ * K + K~ , or J/ip — > K + K~ 
decay. D + — > K S K + candidates are also removed. We use the same rf selection criteria as 
those used in our previously published analyses 0, 0| ■ 

For reconstructed B —>■ fcp candidates, we identify B meson decays using the en- 
ergy difference AE = E^ — E^ m and the beam-energy constrained mass My, c = 



(E™^ m ) 2 — (p c ™) 2 , where E^ m is the beam energy in the cms, and E^ s and p c ™ are 
the cms energy and momentum of the reconstructed B candidate, respectively. The B me- 
son signal region is defined as \AE\ < 0.06 GeV for B° — > (pK^, \AE\ < 0.04 GeV for 
B° -> K + K~K° S , \AE\ < 0.06 GeV for B° -> jf{-^ pj)K^, or -0.10 GeV < AE < 0.08 
GeV for B° -> ^'(-> n + n-r])K° s , and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 for all decays. In 
order to suppress background from the e + e~ — > mm, dd, ss, or cc continuum, we form signal 
and background likelihood functions, £3 an d ^bg^ from a set of variables that characterize 
the event topology, and impose thresholds on the likelihood ratio £s/(£>s + £bg) [IH- The 
threshold value depends both on the decay mode and on the flavor-tagging quality. 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed B° — > fcp decay 0|. We use two 
parameters, q and r, to represent the tagging information. The first, q, is already defined 
in Eq. ((TJ). The parameter r is an event-by-event, MC-determined flavor-tagging dilution 
factor that ranges from r = for no flavor discrimination to r = 1 for unambiguous flavor 
assignment. It is used only to sort data into six r intervals. The wrong tag fractions for the 
six r intervals, wi (I = 1, 6), and differences between B° and B° decays, Awi, are determined 
from the data; we use the same values that were used for the sin 20! measurement [l5j | . 

The decay vertices of B° mesons are reconstructed using tracks that have enough SVD 
hits. The vertex position for the fcp decay is reconstructed using charged tracks excluding 
pions from the K$ decays. The / tag vertex determination remains unchanged from the 
previous publication and is described in detail elsewhere [lif . 

After flavor tagging and vertex reconstruction, we obtain the numbers of B° — > fcp 
candidates, N ev , listed in TablelU Figured] shows the Mb c distributions for the reconstructed 
B candidates that have AE values within the signal region. We use events outside the signal 
region as well as a large MC sample to study the background components. The dominant 
background comes from continuum events. In addition, according to MC simulation, there 
is a small (~ 3%) contamination from BB background events in B° — ► rfKg (7/ — ► p 0r y). 
The contributions from BB events are smaller for other decay modes. The contamination of 
K + K~K$ events in the <pK^ sample (7.2 ± 1.7%) is also small. Finally, backgrounds from 
the B° f (980)Kg decay, which has the opposite CP eigenvalue to (j)Kg, are found to be 
small The influence of these backgrounds is treated as a source of systematic 
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TABLE I: The numbers of reconstructed B° —> fcp candidates used for S and A determination, 
iVev, and the estimated signal purity in the AE-M^ C signal region for each fcp mode. 



Mode 


£/ 


AW 


Purity 


<t>K« 

K + K~Kg 

n'K° s 


-l 

+1(100%) 

-1 


106 
361 
421 


0.64 ±0.10 
0.55 ±0.05 
0.58 ±0.05 



30 
20 [ 
10 \ 



L a) B° ^ ^K° s 



O 
> 

CD n 

c<3 100 h 
o 
o 

d 50 k 

i o 



LU 



b) B° -> K + K"K° S 

+ ++ + ii ++ ^ i ^ + ^ t_±_^ 



150 L c) B°^ti'K° s 
100 - 



50 ; + f i ++ 4 J i- + ++ i + i + i+i + + i + i i - 




5.22 5.24 5.26 5.28 

M bc (GeV/c 2 ) 



5.3 



FIG. 1: The M bc distributions for (a) B° -» (f)K% (b) B° -> K+K~K° S , and (c) £ c 



within the A.E signal region. Solid curves show the fit to signal plus background distributions, and 
dashed curves show the background contributions. The fit functions are the same as those used in 
the previous publication 



uncertainty. 

We determine S and A for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for 
the signal distribution, "P sig (At; S, A, q, Wi, Awi), is given by Eq. ((TJ incorporating the effect 
of incorrect flavor assignment. The distribution is convolved with the proper-time interval 
resolution function _R sig (At) [l5|, which takes into account the finite vertex resolution. We 
determine the following likelihood value for each event: 



Pi 



(1 " /ol) / 

J — ( 



/si g Psig(At') J Rsig(At i - At') 



+ (1 - / S ig)nkg(At')^bkg(A^ - At') 
+ UP ol (AU) 



d(At') 



(2) 
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TABLE II: Results of the fits to the At distributions. The first errors are statistical and the second 
errors are systematic. The third error for the K + K~ Kg mode arises from the uncertainty in the 
fraction of the CP-odd component. 



Mode 


-ZfS (= sin 20i in the SM) 


A (= in the SM) 


<PK S 


-0.96 ± 0.50l#ii 


-0.15 + 0.29 + 0.07 


K + K-K s 


+0.51 ± 0.26 ± 0.05t[J;oo 


-0.17 + 0.16 + 0.04 


n>K s ' 


+0.43 ± 0.27 ± 0.05 


-0.01 + 0.16 + 0.04 



CD 

E 
E 

>s 

CO 

< 

en 
cr 



1 ■ 

0.5 


-0.5 
-1 
1 ■ 

0.5 ■ 
' 

-0.5 ■ 
-1 



B°^K° S 



0.0 < r< 0.5 




0.5 < r< 1.0 



-7.5 -5 -2.5 2.5 
At (ps) 



5 7.5 




-7.5 -5 



-2.5 2.5 
At ( P S) 



5 -2.5 2. 
At ( P S) 



5 5 7.5 



FIG. 2: (a) The asymmetry, A, in each At bin for i? c 



<j>K% with < 



< 0.5, (b) with 



V'K° S 



0.5 < r < 1.0, (c) for B° -» K + K~K% with < r < 0.5, (d) with 0.5 < r < 1.0, (e) for B° 
with < r < 0.5, and (f) with 0.5 < r < 1.0, respectively. The solid curves show the result of 
the unbinned maximum-likelihood fit. The dashed curves show the SM expectation with sin2c/>i = 
+0.731 and A = 0. 



where P Q \(At) is a broad Gaussian function that represents an outlier component with a 
small fraction f Q \. The signal probability / s ; g depends on the r region and is calculated on 
an event-by-event basis as a function of AE and M^ c 0. Pbkg(Ai) is a PDF for background 
events, which is modeled as a sum of exponential and prompt components, and is convolved 
with a sum of two Gaussians i?bkg- All parameters in Pbkg(At) and -Rhk ? are determined by 
the fit to the At distribution of a background-enhanced control sample [17j|; i.e. events away 
from the AE-M^ C signal region. We fix r^o and Am^ at their world-average values (l8j |. The 
only free parameters in the final fit are S and A, which are determined by maximizing the 
likelihood function L = \[ i P^Ati] S, A) where the product is over all events. Table |LT1 gives 
the fit values of — £/<S and A. These results are consistent with the previous results |9( and 
supersede them. We obtain values consistent with the present world average of sin 2<fti = 
+0.731 ± 0.056 in the B° — > K + K~ K Q S and if Kg decays, while a negative value is 
observed in B° — > <pK s decay. 

We define the raw asymmetry in each At bin by A = (A^^ /= _i — Nq£.=+i) / (N q £.=-i + 
N q £ f=+ i), where A r g ^ /=+1 (_ 1 ) is the number of observed candidates with q£f = +1(— 1). 
Figures I2fa-f ) show the raw asymmetries in two regions of the flavor-tagging parameter r. 
While the numbers of events in the two regions are similar, the effective tagging efficiency is 
much larger and the background dilution is smaller in the region 0.5 < r < 1.0. The observed 
CP asymmetry for E>° <pKg in the region 0.5 < r < 1.0 [Fig. [2(b)] indicates the difference 
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from the SM expectation (dashed curve). Note that these projections onto the At axis do 
not take into account event-by-event information (such as the signal fraction, the wrong tag 
fraction and the vertex resolution), which is used in the unbinned maximum-likelihood fit. 

Fits to the same samples with A fixed at zero yield — £fS = —0.99 ± 0.50(stat) for 
B° -> <f)K s , -Z f S = +0.54±0.24(stat) for B° -> K + K~K%, and -£/<S = +0.43±0.27(stat) 
for B° — > rj'Kg. Applying the same fit procedure, we also obtain S = — 0.09±0.26(stat), A = 
+0.18±0.20(stat) for B+ -> <pK + decay and S = +0.10±0.14(stat), A = -0.04±0.09(stat) 
for B + — > 7]'K + decay. Both results for the S term are consistent with no CP asymmetry, 
as expected. 

The dominant sources of systematic error for the B° —>■ <pK s mode are a possible fit bias 
for the input S value near the physical boundary ( j^oo f° r •> the uncertainties in the B° — ► 
K + K~K S and f (9S0)K° s background fractions for S and ±0.04 for .A), in the other 

background fractions (±0.05 for S and ±0.04 for A), and in the vertex reconstruction (±0.02 
for S and ±0.05 for A). Other contributions come from uncertainties in the background At 
distribution, wrong tag fractions, t b o, and Am^. We add each contribution in quadrature 
to obtain the total systematic uncertainty. Systematic uncertainties from these sources are 
also examined for the other modes. We find that the dominant sources are uncertainties 
from the background fractions and from the vertex reconstruction. 

We use the Feldman- Cousins frequentist approach [19J to determine the statistical signif- 
icance of the observed deviation from the SM expectation in B° — > 4>K S . The procedure is 
described in detail elsewhere |2(|. From 1-dimensional confidence intervals for S with A set 
at zero, the case with S = ±0.731 is ruled out at 99.95% confidence level, equivalent to 3.5<r 
significance for Gaussian errors. As a cross check, we varied the selection criteria for the 
B° — > <f)K s decay and repeated the analysis. We find no sizable change in the significance. 

In summary, we have performed improved measurements of CP- violation parameters for 
B° — > <f)K s , K + K~~K S and r)'K s decays. These charmless decays are sensitive to possible 
new CP-violating phases. Our results for B° — > rj'K s and K + K~K S are consistent with 
those obtained for B° — > J /ipK s and other decays governed by the b — > ccs transition. On 
the other hand, a 3.5a deviation is observed for B° — > <pK s . The result suggests that there 
is a large CP-violating phase in its decay amplitude, which cannot be explained by the SM. 
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